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The exciton-polariton propagation in resonant hybrid (isotropic/anisotropic) periodic stacks, with 
misaligned in plane anisotropy and Bragg photon frequency in resonance with Wannier exciton of 
2D quantum wells, is studied by self-consistent theory and in the effective mass approximation. 
The optical tailoring of this new class of resonant Bragg reflectors, where the structural periodicity 
of a multi-layer drives the periodicity of the in-plane optical C axis orientation, is computed for 
symmetric and non-symmetric elementary cell by conserving strong radiation-matter coupling and 
photonic band gap. We will demonstrate, by selected numerical examples, that the behavior of the 
so called intermediate dispersion curves (IDC), that drop between upper and lower branches of the 
lowest energy band gap, are strongly dependent from in-plane C axis orientation. Therefore, we 
guess that this class of hybrid meta-materials is promising for new trapping light optical devices 
based on IDC behavior. 

PACS numbers: 78.67.De, 71.36.+c,42.70.Qs,78.66.Fd 


I. INTRODUCTION 

The manipulation of the optical properties in order to 
obtain coherent radiative coupling among a collection of 
emissive species was proposed by R.H. Dickai in 1954. 
The optical behavior of the composite exciton-polariton, 
propagating in a resonant ID Bragg quantum wells, was 
studied by E.L. Ivechenko et all.-' in 1994, while the sat¬ 
uration effect, that transforms the super-radiant mode 
into a polaritonic Bloch band, usually obtained for rather 
large quantum well number (N > 200), was theoretically 
investigated by K. Cho et all—in 2002 and measured 
by J.P. Prineas at all.— in the same year. 

It is well known that Wannier exciton is intrinsically 
a giant dipole excitation, since the microscopic dipole 
transitions of the elementary cells of a semiconductor are 
coherently coupled in a volume determined by its Bohr 
radius^. Moreover, the coupling between the center-of- 
mass motion, of the exciton, and the photon wave vec¬ 
tor, that introduces the optical spatial dispersion proper¬ 
ties in the polariton model^, in principle can give strong 
radiation-matter coupling also in a single ’’high quality” 
quantum well, if the polariton splitting energy is greater 
than the total broadening of the system. 

When a cluster of N emitters is arranged periodically 
in an homogeneous dielectric background, with emission 
wavelength equal two times the spatial periodicity, a new 
collective coherent state (super-radiance mode) is present 
in the system, also called resonant photonic Bragg reflec¬ 
tors (RPBR). Since the strength of N oscillators concen¬ 
trate into the super-radiant mode, the radiation-matter 
interaction increases by increasing the number of quan¬ 
tum well, and therefore polariton splitting energy can 
overcome the total broadening of the system (the so 
called weak/strong coupling transition^dil). 

Recently, the RPBR becomes an important tool for 
the control and manipulation of light. In fact, the stop¬ 
ping, storing and releasing light in RPBR was recently 
proposed by Yang et all— the former effect is based 


on the parametric manipulation of photonic band struc¬ 
ture, and practical limitations for its realization, due to 
the large number of quantum wells present in the clus¬ 
ter (N > 200), was also discussed in the same reference. 
Moreover, a multilayer system, where exciton energy is 
close to the transmission peak at high energy side of the 
stop band was studied experimentally by Askitopoulos 
et all— i as a new tools for tailoring of light-matter in¬ 
teraction, and promising for polariton lasing, which re¬ 
quires no population inversio n 12 : 13 . At variance of micro¬ 
cavities, where the number of quantum wells are strongly 
limited, the RPBR, obtained by alternating isotropic ma¬ 
terials with different background refraction indices and 
quantum well resonances, close to the Bragg frequency, 
require a number of elementary cells (N=30-60), greater 
than that present in the BDR microcavity, but no so 
large as those required in order to build up the former 
polaritonic gap^£. 

The study of anisotropic photonic crystals, obtained 
by alternating different transparent uni-axial layers (or 
isotropic/anisotropic layers), is largely present in the re¬ 
cent literature due to the new interesting optical proper¬ 
ties, namely: negative refraction—it— , omni-directional 
transmission^ and giant transmission effects*— Re¬ 
cently, in a periodic system obtained by alternating uni¬ 
axial bilayers, with different orientation of the in-plane 
optical C axis, it was demonstrated that band gap can be 
opened inside the Brillouin zone, where band curvature 
indicates the presence of negative refraction— 

The present work is devoted to study the optical re¬ 
sponse of a new class of resonant hybrid Bragg reflectors 
(RHBR), where the periodicities, due to the background 
dielectric constant modulation, and different orientation 
of the optical C axis on the xy plane are present, and 
that will be able to combine the optical properties of the 
isotropic resonant quantum confined systems, namely: 
super-radiance, strong radiation-matter interaction, high 
non-linearity, with the interesting optical properties ob¬ 
served in anisotropic multilayers, in particular: high pho- 
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elementary cell (N=l) 


FIG. 1: (color online) ID hybrid (isotropic/anisotropic) 
multilayer. Black (brown online) layers are the quantum 
wells, light grey (yellow online) layers are the barriers of 
the quantum wells and dark grey (orange online) layers 
are the uniaxial bilayers. The dashed arrow indicates 
that C axis could not be parallel to the x axis. 



FIG. 2: (color online) Full width at half maximum of 
the reflection spectra of a resonant Bragg quantum 
wells at normal incidence in an isotropic bulk as a 
function of the number N of elementary cells. 
Parameter values are given in Sec. II of the paper. 


ton density of states at absolute photonic gap, its build¬ 
ing up for a rather moderate number of elementary cells, 
and in general stronger flexibility in the optical tailoring. 
Finally, optical properties of RHBR will be compared 
with those shown by analogous RPBR. 

Notice, that the recent improvement on plasma- 
assisted molecular-beam epitaxyi^ have allowed to study 
multilayer with in-plane optical C axis, and, even if the 
goal of the total control on C axis orientation in a sin¬ 
gle multilayer is not very close, the present study of the 
optical response in RHBR should be a stimulation for 
studying the optical properties of systems, interesting for 
fundamental and applications, as: i) meta-material topo¬ 
logical insulator—, ii) repulsive Casimir forced and iii) 
the so called chiral meso-molecules— 

The aim of the present work is twofold: a) to inves¬ 
tigate, for selected physical parameter values, exciton- 
polariton propagation in a RHBR stacks of N-Bragg 
quantum wells under strong radiation-matter interaction 
for symmetric and asymmetric elementary cell, and, b) 
to study the dispersion curves of the former systems, in 
order to investigate the behavior of the so called interme¬ 
diate exciton-polariton dispersion curves (IDC), in the 
lowest energy photonic stop band. The latter point is 
particularly interesting for light storing as shown in Ref. 
10 for isotropic RPBR. The exciton-polariton propaga¬ 
tion is computed in a finite RHBR, composed by N el¬ 
ementary cells, where a 2D quantum well is located in 
an isotropic A/4 layer, and the other layer is composed 
by uniaxial bilayers, with in-plane C a axis (a = L,R), 
located at left and right side of the uniaxial bilayer, as 
schematically reported in the Fig. 1 for the case N=1 el¬ 
ementary cell. Notice, that in all the present calculation 
we consider a cluster of N elementary cells at A/2 for a 
total of N+2 quantum wells, due to the insertion of two 
surface quantum wells (see Fig. 1). 


The study of photonic dielectric gap and super¬ 
radiant exciton-polariton mode building up, as a func¬ 
tion of quantum well number N in a RHBR, is per¬ 
formed by self-consistent calculation and in effective mass 
approximation^. Notice, that in the present study the 
use of non-local self-consistent framework is mandatory 
for a correct estimation of the radiative self-energy, whose 
imaginary part, complemented with non-radiative ho¬ 
mogeneous broadening (usually derived from the experi¬ 
ments), allows to asses the reaching of the strong cou¬ 
pling regime in the radiation-matter interaction^ and 
gives a quantitative estimation of the absorbance spec¬ 
tra in the so called ’’high quality” quantum wells, where 
the non-homogeneous broadening should be negligible. 
In order to obtain an heavy-hole Wannier exciton opti¬ 
cally isotropic, in an otherwise anisotropic super-lattice, 
we have to minimize the envelope function penetration 
depths into the uniaxial layers and also the effect of the 
image potential on the exciton dipole, due to the presence 
of different uniaxial layers at both sides of the quantum 
well (see Fig. 1), as discussed by V. Trongin et all—. 
The former requirement is accomplished by considering 
an heavy-hole Wannier exciton perfectly confined into 
a 2D quantum well, and clad between two rather large 
(with respect to the exciton Bohr radius) isotropic bar¬ 
riers, while the latter requirement is obtained by tak¬ 
ing the background dielectric constant equal for well and 
isotropic barriers (£b = £ w ), and moreover not very differ¬ 
ent in value from the orthogonal (e±) and parallel (ey) 
dielectric constants of the uniaxial layers. Now, since 
the main approximations, embodied in the model calcu¬ 
lation, are briefly discussed, let us summarize the plan 
of the work. In Sec. II the procedure for computing 
the optical response of a RHBR, for symmetric and non- 
symmetric elementary cell, is presented, and also the be¬ 
havior of the exciton-polariton propagation in RPBR, for 
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homogeneous and periodically modulated dielectric back¬ 
ground, are briefly summarized^. In Section III the op¬ 
tical response for a finite RHBR with in-plane C axis is 
computed for symmetric and non-symmetric elementary 
cell. The computation is performed for N=32 elementary 
cells (and N=34 quantum wells) where strong coupling 
regime and photonic band gaps are observed for selected 
values of physical parameters. Their absorbance proper¬ 
ties and the behavior of radiative homogeneous broaden¬ 
ing and non-radiative broadening, particular important 
in optical applications, will be discussed. In Section IV 
the dispersion curves of former systems will be computed 
(N —> oo) for non-radiative broadening value Tnr —> 0. 
The behavior of the IDC in resonance and out of reso¬ 
nance with Bragg frequency will be discussed at variance 
of the correspondent isotropic casei The conclusions will 
be summarized in the Sec. V. 


II. THEORY 

The study of photonic dielectric gap building up and 
super-radiant exciton-polariton propagation, as a func¬ 
tion of quantum well number N, is performed by self- 
consistent calculation, and in the effective mass approx¬ 
imation along the line of Refs. 4,6 and 7. A standard 
method based on Green functions and transfer matrix 
approach is adopted for numerical calculation, and all 
the formula are explicitly given in the appendix. Now, 
let us shortly remind the radiation-matter properties of 
N Bragg quantum wells in an isotropic background, com¬ 
puted for the same physical parameter values that will be 
used for the isotropic slab of the hybrid photonic crystal 
of Sec. III. 

In the present calculation we have used quan¬ 
tum well parameter values, close to those of the In- 
GaAs/AlGaAs/GaAs(001) semiconductor materials, as 
for istance: total mass M = 0.524, Bohr radius qr = 
8.047nm and background dielectric constants: e w = Eb = 
10.24. In all the cases the used well width is L w = 10 nm. 

The exciton transition energy is: E ex (K\\) = E ex ( 0) + 
h 2 K^/2M, where Ky is the in-plane center-of-mass wave 
vector and, in the present calculations, Ky = q x while 
E e: „(0) = 1.418eU is the variational energy of the exciton 
confined between infinite potential barriers. The homo¬ 
geneous non-radiative broadening T nr is chosen close to 
that is experimentally observed in high quality quantum 
wells at rather low temperature: Tnr = 0.25meU and 
T « 10 K. 

Finally, in order to minimizes the Fabry-Perot oscil¬ 
lations, due to the vacuum/semiconductor background 
dielectric mismatch, that strongly perturb exciton- 
polariton propagation, the ID cluster is clad between 
two isotropic semi-infinite bulk materials with Eb back¬ 
ground dielectric constant. Notice that this effect is usu¬ 
ally obtained experimentally by an anti-reflection coat¬ 
ing deposited on both the surfaces of the sample. Now, 
let us briefly summarize the different regimes of exciton- 


polariton propagation in an ID cluster of N quantum 
wells at A/2 separation embedded in a homogeneous bulk 
background (Z > 0) with refraction index Eb = 10.24. 

In Fig. 2 the full width at half maximum of the reflec¬ 
tivity is shown as a function of the number N of quantum 
wells, for resonant wavelength: A = 2i tc/ (u} ex y/Eb) = 2d, 
where d is the periodicity and hui ex = E ex ( 0). Three 
different regimes can be observed^, namely: the linear 
behavior, due to the super-radiant mode, at low N quan¬ 
tum well number (N = 1 -y 80), the saturation zone, 
obtained for large quantum well number (N > 130), 
where multi-modes polariton propagation are observed 
(that for N —> oo will merge in forward and backward 
polariton Bloch propagation) and an intermediate zone 
(80 < N < 130), where the system changes from super¬ 
radiance to multi modes behavior. Notice, that in the 
zone of super-radiant mode the system can change its 
behavior from weak to strong radiation-matter coupling 
when the polaritonic splitting energy, that is a direct 
function of the number N of quantum wells present in 
the stack, becomes greater than the total (radiative plus 
non-raditive) homogeneous broadening (Tj = T nr+Tr). 

At variance of the behavior shown in Fig. 2, where 
the polaritonic stop band is due to the Bragg periodic¬ 
ity of the dispersive component of the exciton suscep¬ 
tibility (N > 130), the photonic stop band due to the 
isotropic dielectric background modulation, shows a more 
fast convergence (for N ze 30 -y 60) and broader stop 
band in energy also for rather low dielectric contrast 
Ae = e 2 — Ei. In fact for £2 = £b and £1 = 8.24, the 
stop band width of an isotropic Bragg reflector in reso¬ 
nance with hui Q = E ex { 0), computed by the equation: 
AA = (4/7r) A 0 sm _1 (n 2 - ni) / (n 2 + ni) (jn = ^/eI), 
is Aw 0 « 81meF , while for homogeneous background 
we obtain: Aw n ss 19meV\ In conclusion, while the 
exciton-polariton propagation in a resonant Bragg sys¬ 
tem, with homogeneous dielectric constant background, 
shows three different well separated regimes as a func¬ 
tion of N (as shown in Fig. 2 for N = 1 -y 250), the 
exciton-polariton propagation in a resonant Bragg sys¬ 
tem, with background dielectric modulation, show a zone 
of N-values (N = 30 -y 60) where super-radiant and pho¬ 
tonic gap regimes are both present in the same cluster. 


III. OPTICAL RESPONSE OF PHOTONIC 
HYBRID CLUSTER 

In the present calculation we consider the optical re¬ 
sponse in a RHBR, composed by N symmetric elementary 
cells with z axis along the periodicity and (x,y) in-plane 
coordinates, for incident wave polarized S or P, and in¬ 
plane optical C axis of the uniaxial slab («l = cur = a) 
as shown in Fig. 1 for N = 1. 

As discussed in the introduction, the uniaxial refrac¬ 
tion index values are chosen very close to the isotropic 
refraction index rib, namely: n± — rib = 3.2 and n\\ = 3.5 
respectively. The former choice emphasizes the role of 
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FIG. 3: (color online) P-polarized optical response, at normal incidence, of a RHBR (N=32), computed for different 
a values. Reflectivity: dashed line (red online); transmission: dotted line (bleu online) and absorbance: solid line 

(green online). 



energy (eV) 

FIG. 4: (color online) Transmission x-component ( T pp ): 
dashed line (red online); y-component ( T ps ): dotted line 
(bleu online) and reflection x-component (I? pp ): dark 
solid line (green online); y-component (R ps ) light solid 
line (grey online), computed for the system of Fig. 3c. 


the orientation of the optical C axis on the xy plane, 
and makes simpler to compare the total optical response 
in RHBR with that observed in equivalent RPBR sys¬ 
tem. Notice that the tuning between the exciton en¬ 
ergy and the Bragg energy of the photonic stop band 
in uniaxial periodic multi-layers with symmetric cells 
can be computed straightforward for a = n 7r/2-values 
(n = 0,±l,±2,...), since the dielectric tensor of uniaxial 
slab becomes diagonal, and therefore, ordinary and ex¬ 
traordinary wave components do not interact. For a clus¬ 
ter of N=32 symmetric elementary cells, strong coupling 
regime is observed^, and in order to obtain the Bragg 
energy hui 0 , of the reflection stop band at normal inci¬ 
dence, in resonance with exciton energy (hco 0 ~ .E^O) = 
1.418614), for incident P polarized wave we will have 
to compute the A/4 isotropic thickness for wavelength: 


A = 2nc/ ( u> ex rib ), while the total thickness of theA/4 
uniaxial layer is given for wavelength: A = 2irc/ ( U) ex nu ). 
Notice that, due to the optical symmetry of the system^, 
the incident P-polarized optical response, with a ori¬ 
entation of the C axis, shows the same optical spectra 
than that computed for S polarized incident wave and 
a' = (7t/2) — a orientation. 

In Fig. 3a the total optical response at the normal 
incidence for P polarized incident wave is shown. We 
observe a reflection central deep, very close to the bare 
exciton energy, in resonance with a rather intense ab¬ 
sorbance peak (Ia ~ 88%), and a rather broad photonic 
stop band of A u)(q x = 0) = lOlmeP (full width at half 
maximum)^d£. Moreover, since the dielectric modula¬ 
tion of the multilayer is due to the £|| dielectric constant 
and to the dispersive component of the exciton suscepti¬ 
bility, while for the ordinary wave the multilayer results 
homogeneous (ej_ = £&), the band gap disappears in the 
spectrum computed for C axis orientation along the y 
axis (a = 7 t/ 2) as is shown in Fig. 3b. In fact, only 
the normal exciton-polariton optical spectra are present, 
where the absorbance shows a double peak with two max¬ 
ima separated in energy of about 1.7meV. Notice, that in 
the present kind of RHBR the reflection intensity can be 
tuned from 0% to 100% by changing the orientation of 
the C axis from a = 0 to a = 7r/2, while the line-shape 
remain unchanged and the transmission spectrum, close 
to the border of the reflection stop band, maintains its 
high intensity value (It ~ 100%). In fact, the optical 
response for P incident polarization and a = 7t/ 4 , where 
an equal contribution from ordinary and extraordinary 
waves is expected, is shown in Fig. 3c. The main optical 
feature of Fig. 3c is the reflection line-shape that is very 
similar to the stop band shape except for the intensity 
that now is close to 50% of reflection, instead of about 
100% observed in Fig. 3a. Analogously, for the trans¬ 
mission line shape in the gap energies, that shows a com¬ 
plementary behavior with respect to the reflection one, 
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(c) P-polarization aj, = 45°, aj? = 0° 


(d) S-polarization ax = 45°, or = 0° 


FIG. 5: (color online) Linear-polarized optical response, at normal incidence, of a RHBR (N=32), computed for 
different ax values. Reflectivity: dashed line (red online); transmission: dotted line (bleu online) and absorbance: 

solid line (green online). 


with the 50% of the transmission intensity. Notice that, 
the former property is, almost verified also by removing 
the condition e± = Eb (for instance^): ex = 8.24. 

In order to go a bit deeper in understanding the optical 
behavior of the former system, let us analyze the optical 
response as a function of its polarized components. In 
Fig. 4 for P polarized incident wave the reflection and 
transmission in-plane components are shown. While the 
total reflection spectra, for photon energies in the stop 
band, are obtained by a rather equal contribution of R pp 
and R ps components, the contribution of the T pp compo¬ 
nent at total transmission spectra is essentially localized 
at low energy side of the stop band, while the T ps contri¬ 
bution is at high energy side. Obviously, for S polarized 
incident wave the pp and sp transmission components 
must be exchanged. 

In conclusion, in symmetric RHBR, with N = 32, the 
absorbance spectrum clearly show two components due 


to IDC bands, and moreover, for selected values of the 
physical parameters, it is possible to determine the orien¬ 
tation of the in-plane C axis by reflection and/or trans¬ 
mission intensity measurements, and also the polariza¬ 
tion of the incident wave by analyzing the polarization 
of the transmission peaks at low and high energy side of 
the photonic stop band. Usually, the optical tailoring of 
RHBR require the solution of the so called inverse prob¬ 
lem, while in the present case these informations can be 
deduced directly from a quantitative analysis of the op¬ 
tical response. Now, let us consider a bit more complex 
system than before obtained by substituting the uniaxial 
layer with an uniaxial bilayer of equal thickness, with the 
optical Cp axis (with /3 = L,R) oriented along x-axis for 
the right layer (an = 0) and in-plane for the left layer 
(0 < ax < 7 t/ 2) respectively. Notice, that in this case 
the elementary cell is asymmetrical, except in the limit 
of ax —t 0, and therefore we will have to distinguish be- 
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tween forward and backward optical response. 

In the present computation a different choice than be¬ 
fore is adopted for e± in order to accomplish the follow¬ 
ing two conditions: e± < £b = £ w < £y where ex = 8.24 
and (ey + £j _)/2 ss e 6 , moreover an equal A/4 thickness 
value, computed for background dielectric constant ej,, is 
chosen for isotropic and anisotropic layers, in resonance 
with bare exciton energy (fio> 0 (0) ~ E ex ( 0)). First of all, 
let us to compute total reflectivity and adsorption for S 
and P incident waves for a multilayer with N=32 ele¬ 
mentary cells in the limit of symmetric elementary cell 
(ax = au —> 0). In this case the system is similar to that 
discussed before, except for the thickness of the uniaxial 
bilayers that now is in resonance neither with A/4 com¬ 
puted by n\\ nor with that computed by n± values. The 
numerical results are shown in Fig. 5a and 5b for P and S 
incident polarization respectively, where two well shaped 
reflection stop bands can be observed at opposite energy 
sides of the exciton absorption peak, and moreover in 
the absorption spectra a double peak is present. Now, 
by increasing the ax value (with ctR taken constant and 
equal zero) the forward optical response shows that the 
two stop bands move in the opposite directions in energy, 
towards the exciton energy value, that will be reached 
for ax = 7r/4 as shown in Fig. 5c and 5d for P and S 
incident wave polarization respectively. Notice that for 
the former ax value a complete mixing between ordinary 
and extraordinary waves are expected, and a well formed 
photonic stop band (Aw(0) ~ 57me V) in resonance with 
bare exciton energy are present for both the polariza¬ 
tions. Notice, that at first sight the reflection spectrum 
of Fig. 5c could be interpreted as a very broad reflection 
band (Aw(0) ~ 114me?;) where the Bragg energy seems 
not in resonance with the exciton energy, but that this 
interpretation is not correct can to be demonstrated from 
the location in energy of the two side bands^ present in 
the exciton absorbance spectrum (not reported here) and 
from the dispersion curves that will be discussed in the 
next section. 

Finally, for ax = 7 t/ 2 two photonic band gaps, with 
Bragg energies in resonance with the exciton energy 
(Aw(0) ss 103meu), are shown in Fig. 6a and 6b for 
both the polarizations. Notice, that in Figs. 6a and 
6b the two spectra show exactly the same line shapes, 
due to the ax = 7t/2 condition^, except for photon en¬ 
ergies very close to the exciton transition energy, where 
exciton absorptions show rather different intensities, but 
analogous line shapes, for S and P incident wave polar¬ 
izations as shown in Fig. 7. Since the forward absorption 
spectra show the same rather small absorbance intensity 
for S and P incident wave polarizations under the con¬ 
dition Tnr —> 0 (not reported here), we can conclude 
that in the present case the homogeneous non-radiative 
broadening affects essentially P polarized incident wave 
spectrum, while gives a very small contribution to S- 
polarized one. The former effect is due to the asym¬ 
metry of the elementary cell, as it could be easily verified 
by observing that S and P absorption spectra exchange 


their intensities, by exchanging forward with backward 
in the exciton-polariton propagation. Moreover, the for¬ 
mer interpretation is in complete agreement also with 
qualitative considerations derived from the computation 
performed in an equivalent isotropic RPBR. 

In conclusion, a quantitative estimation of non- 
radiative homogeneous broadening can be obtained from 
experimental optical response performed in RHBR by 
adopting an asymmetric elementary cell. Notice that, 
the former interesting result is easily generalized to the 
RPBR with asymmetric elementary cell. 

Now, before to discuss the dispersion curves (N —> 
oo) of the former system, for sake of completeness, let 
us show the forward optical response of P incident wave 
polarization of Fig. 6a, reported in a much larger range 
of photon energies where also higher energy stop bands 
are present (see Fig. 8). In this case the absorbance 
peak, correspondent to the IDC, at the Bragg energy of 
the lowest energy stop band, is also clearly shown. For 
S incident wave polarization we should observe the same 
line shape as for P polarization, except for intermediate 
exciton-polariton absorbance peak that in this scale of 
intensities should be hardly observed. 


IV. EXCITON-POLARITON DISPERSION 
CURVES 

In order to go a bit deeper in understanding the 
exciton-polariton propagation in the former hybrid pe¬ 
riodic system with asymmetric elementary cells the com¬ 
puted dispersion curves (for N —> oo and T nr —> 0) are 
shown in Fig. 9a for ax = 7 t/ 4. Notice, that the optical 
response for a cluster of N=32 asymmetric elementary 
cells are reported in Fig. 5c and 5d for incident waves 
polarized P and S respectively. The photonic energy gaps 
drop inside the Brillouin zone, as clearly observed in the 
second and third gap of the picture, and this is a rather 
general property of RHBR with in-plane -axis, due to the 
anisotropic properties of the elementary cell as discussed 
in Ref. 17. In the lowest energy gap of Fig. 9a an IDC, 
close to the unperturbed exciton energy, is also reported. 
This feature, that is present also in isotropic Bragg quan¬ 
tum wells^, characterizes the exciton-polariton propaga¬ 
tion and is promising for optical applications; in fact, its 
rather constant behavior in energy supports zero polari- 
ton group velocity, therefore a trapped light in the RHBR 
could be obtained if this property will be maintained in 
whole the Brillouin zone as discussed, in Ref. 10, for 
isotropic RPBR. 

In Fig. 9b the first gap is shown in an enlarged 
energy scale. Notice, that the first photonic gap 
(Aw(0) ~ 57 meV), is in very good agreement with that 
observed in the cluster of Fig. 5d, for S incident wave po¬ 
larization (dark curves), while for P polarization (light 
curves) a further high reflection zone, at lower photon 
energies, is present. Notice, that the lowest energy dis¬ 
persion curve interacts, with the exciton-polariton IDC, 
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(a) P-polarization tij, = 90° a:= 0° 


(b) S-polarization m = 90° an = 0° 


FIG. 6: (color online) Optical response, at normal incidence, of a RHBR (N=32). Reflectivity: dashed line (red 
online); transmission: dotted line (bleu online) and absorbance: solid line (green online). 



energy (eV) 

FIG. 7: (color online) Absorbance spectra of the system 
of Fig. 6a (dashed line) and Fig. 6b (solid line) in 
enlarged scale of intensity (A = 0 4- 0.01). 


giving rise to the asymmetric behavior of the optical re¬ 
sponse of Fig. 5c, as discussed in the former section. 

Very close to the boundary of the Brillouin zone, the 
lower and upper band of the absolute photonic gap show 
a strongly deformed behavior, with respect to the nor¬ 
mal parabolic one, that should give high polariton den¬ 
sity of states, and this property is related with the so 
called ’’gigantic transmission band edge resonances^”. 
Moreover, two IDC are present in the gap, that are very 
close to the unperturbed exciton energy, except for quasi¬ 
momentum K-values, in corrispondence of the direct gap, 
where the interaction with lower and upper photonic dis¬ 
persion curves removes the degeneracy and therefore they 
show an energy separation as large as Aw ps 5.5 meV. 
The former energy splitting between the two IDC are due 



1 1.5 2 2.5 3 3.5 4 4.5 5 

energy (eV) 

FIG. 8: (color online) The same optical response of 
Figs. 6a shown in a wider energy scale. 


to the repulsion between the high density of states of the 
upper photonic dispersion curve of the gap with the lower 
energy IDC that show the same symmetry, and between 
the lower dispersion curve with the upper energy IDCii. 
Notice, that for a = n/2, where band structures become 
degenerate for P and S wave polarizations^, these two 
curves recombine in one IDC two fold curve (see the in¬ 
set of Fig. 9b), that is rather constant in energy in the 
whole Brillouin zone, since the interaction between upper 
and lower photonic branches and the two IDC becomes 
strongly reduced due to the increasing of the energy gap 
(Aw(0) changes from about 57meV to about 90meV). 
Moreover, the band gap moves to the border of Brillouin 
zone and shows a normal quadratic behavior . There¬ 
fore, the possibility of opening or closing the window of 
Aw ps 5.5meV by different orientation of the left C axis 
is very promising for studying the light storing effect in 




































FIG. 9: (color online) Dispersion curves of a periodic system with the same parameter values of the former Fig. 3C, 
computed in the limits: IV —> oo and Tnr —> 0. Figs, b and c are shown in a reduced energy range around the 
Wannier exciton energy, x-polarized eigenvalues are represented by light (mauve online) dots and y-polarized by 

dark (bleu online) dots. 




FIG. 10: (color online) Dispersion curves of a periodic system, with the same cell geometry of the former Fig. 9a, 
but computed by using £|| = 13 and e± = 2.25. Fig. b is shown in a reduced energy range around the Wannier 
exciton energy, x-polarized eigenvalues are represented by light (mauve online) dots and y-polarized by dark (bleu 

online) dots. 


RHBR. In fact, at variance of an analogous isotropic case 
discussed in the Ref. 10, where the switch on/off is ob¬ 
tained by dynamical Stark effect, in the present system 
the switch is obtained by a linear renormalization of the 
photonic gap, and this is the most interesting result of the 
present section. Moreover, it is well known that dynami¬ 
cal Stark effect decreases the exciton oscillator strength, 
and therefore the system can undergo an unwanted tran¬ 
sition from strong to weak coupling regime. 

In principle, to remove the degeneracy between the 
two IDC on the energy gap should be obtained also 
for non-normal incidence propagation in the RHBR. In 
fact, let us consider the dispersion curves for non-zero 
in-plane wave vector [q x = q\\ ^ 0) computed for inci¬ 


dent angle 8 = 10°. In this case, while the energy 
shift of the exciton center-of-mass is negligible small 
(h 2 q x /2M ss 6.7 fieV), the shift, due to the longer path 
of the light in the multilayer, is clearly observed for 
= 7 t/ 4 (see fig 9c). In fact, the direct polaritonic 
gap changes from Aw(0) ~ 52 meV to A uj{q x ) « 42meR, 
and the top of the lower photonic band shift in energy of: 
fo[u v (qx) — w„(0)] = 33.3 meV. Notice, that in this case 
only three peaks are present in the optical spectra if the 
relationship A ui(q x ) A uj p is verified; moreover, while 

the energy splitting between the upper IDC and lower 
polariton branch interaction is about Auj p « 10 meV the 
lower IDC seems rather unperturbed, since the upper po¬ 
lariton branch is rather higher in energy. Notice, that the 
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former result is in rather good agreement with analogous 
result obtained in isotropic RPBR in the Ref. 11 where 
the interaction is between upper photonic branch and the 
lower IDC. 

Finally, we have computed the dispersion curves for 
the same structural system discussed before, but with 
stronger dielectric contrast in the uniaxial bilayer (ey = 
13 and £j_ = 2.25). In Fig. 10a more pronounced en¬ 
ergy gaps are present inside the Brillouin zone, there¬ 
fore, for exciton energy close to one of these photonic 
gaps an anomalous exciton-polariton propagation could 
be presen t 16 ’ 17 ; moreover, also in this case, due to the 
shift of the photonic band gap, two IDC, close to the top 
of the first optical valence band, are observed. In Fig. 
10b the first energy gap is shown in an enlarged energy 
scale and with this resolution an indirect absolute pho¬ 
tonic gap is observed (A u>(q x = 0) ~ 119.2meV). In a 
further enlarged energy scale (see the inset of Fig. 10b), 
the former pattern looks rather symmetric with respect 
to the maximum of lower photonic band energy. Two 
dispersion curves are present in the photonic gap with 
strong different behaviors, namely: one dispersion curve 
is a rather unperturbed and shows a negligible group ve¬ 
locity till very close to the boundary of the Brillouin 
zone!, while the second curve shows strong distortion 
due to the repulsion with the top of the lower photonic 
band, its group velocity shows different sign with respect 
to the maximum of the dispersion curved. Since also in 
the present case A u)(q x = 0) A uj p , a pattern of three 
polariton dispersion curves is confirmed in the interaction 
zone!!. 

In conclusions, hybrid isotropic/anisotropic resonant 
photonic crystals seems a rather promising class of meta¬ 
materials, with respect to the isotropic multilayers, be¬ 
cause strongly enlarges the possibilities of the optical tai¬ 
loring for fundamental studies and device applications!!. 


V. CONCLUSION 

In the present work, the non-local optical response of 
a RHBR cluster with in-plane C axis, composed by N 
symmetric and asymmetric elementary cells, is computed 
by self-consistent calculations and in the effective mass 
approximation. The non radiative homogeneous broad¬ 
ening value, that in the model calculation is the only 
parameter derived from the experiments, is taken coher¬ 
ently with those adopted for the so called ’’high quality” 
quantum wells. 

The optical response of a RHBR composed by symmet¬ 
ric elementary cells shows strong radiation-matter cou¬ 
pling and, for selected values of the physical parameters, 
allows to extract, from reflection and/or transmission 
intensities, the orientation of the in-plane C axis, and 
from polarized transmission spectra the polarization of 
the incident wave. Notice that all these information can 
also be derived from the solution of the inverse problem. 
Moreover, the problem of tuning exciton transition en¬ 


ergy and Bragg energy of the stop band in a cluster of 
hybrid multi-layers is also addressed. 

The RHBR with asymmetric elementary cell are ob¬ 
tained by substituting the uniaxial layer of the former 
system with uniaxial bilayer with different orientation of 
the optical axis (see Fig. 1). We show theoretically that, 
for special values of in-plane C axis orientations, the ab¬ 
sorption spectra of this system can give a direct quan¬ 
titative estimation of non-radiative homogeneous broad¬ 
ening of the exciton-polariton in ” high quality” quantum 
wells. Notice, that the former parameter, as underlined 
before, strongly affects the results of the self-consistent 
calculation of the optical response in RHBR. 

Finally, the computation of the dispersion curves in 
RHBR, with asymmetric elementary cell, has allowed to 
investigate the exciton-polariton propagation and/or lo¬ 
calization, due to the presence of two IDC in the lowest 
energy gap, very close in energy to the unperturbed ex¬ 
citon dispersion curve, that show a strong different be¬ 
havior as a function of C axis orientation. We guess that 
the former property is promising for optical device real¬ 
izations. 


Appendix A: EXCITON-POLARITON MODEL IN 
A 2D QUANTUM WELL 

The study of exciton polariton propagation, in semi¬ 
conductor materials, requires the solution of the Maxwell 
equation: 

V x V x E(i?;w) = q 2 £bE(R;uj) + 47tP ex (R;oj) (Al) 

where q = w/ c, P ex (R;oj) is the nonlocal linear polariza¬ 
tion vector: 

L w /2 

P ex {R-uj) = q 2 S ex ^(R) f dR'y*{R')-E{R'-,u) 

-i„/2 

(A2) 

and ’F is the exciton envelope function. 

In isotropic material we can choose Cartesian coordi¬ 
nates with the Z axis perpendicular to the surface plane 
(X,Y) and, for the cylindrical symmetry of the quantum 
well, we can adopt mixed coordinates (K\\,R) for the 
motion of the exciton centre-of-mass. In this case the 
exciton function can be written as: 

^(r,R) = Ntp(r)e iK " R " (A3) 


where r = \Jp 2 + \z e — znf ■ For a 2D heavy hole Wan- 
nier exciton, perfectly confined between infinite potential 
barriers, <p(r) ca be written as a two subbands variational 
envelope function!!: 

<p(r) = N ex cos cos e ~ p/aex (A4) 
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As it is well known, the solutions Ep (/? = x, y) of Eq. 
(Al) can be obtained by the solutions Ep of the corre¬ 
sponding homogeneous part {R ex = 0) combined with a 
solution of the heterogeneous equation (see, for instance, 
Refs. 4 and 6), namely: 

Ep{u,Z) = E°p(uj,Z) - q 2 pS ex (u,q x ) 

L w /2 

x f dZ"GW)**(Z") 

— L w /2 
L w /2 

x f dZ' 'H{Z')Ep(cu : Z') (A5) 

-L m /2 

where Gw^ is the Green function. Notice that, we have 
embodied all the different normalization constants^, of 
the Green function, into the coefficient g 2 that assumes 
the values: q x = k^/£ w and q 2 = w 2 /c 2 respectively. 
Therefore, the Green functions become the same for both 
x and y components: 


G w (z,z') = 


i2q u 


(A6) 


All the energy dependence is contained in the quantity: 

So(co) 


Sex ) — 


E 2 (K ||) — h 2 ui 2 — i2huiT nr 


(A7) 


where S 0 (uj) = 47 TghE K e 2 /wm 0 and E (K ||) the total 
exciton energy peak: 


E {K\\) — E gap E ex + 2M^\\' (A8) 


Ek is the Kane energy ( Ek= 23eV in GaAs based 
semiconductors), m 0 the elctron mass and g the spin- 
degeneracy that, in the present calculation, it has been 
taken as a constant value (g = 1). 

Finally, by solving the Lippmann-Schwinger 
equatior^^, we can obtain the electric field into 
explicit form (with the factor exp [iq x X] suppressed): 


Ep( cc; Z) = [e ik - z - g w (Z) & f>(w; q x ) q> ex {q x ) 
+B [P) \ e ~^ - g w ( Z ) q x ) Vex(-qx) 1 (A9) 


ip ex is the Fourier transform of exciton envelope function: 


L w /2 

q> ex (±k w ) = f 4/(r = 0, Z) e ±iqxZ dZ (A10) 

-L w /2 


tyex {Eq x ) 


Ae/j: k fx . / k w L w \ 

kw [kl x - kl) Sm ) 


(All) 


We define now the integral function g w (Z) computed 
at the slab interfaces (Z = ±L W / 2): 


L w /2 


g w (Z) = J dZ"G w (Z,Z")** ex (r = 0,Z") (A12) 


— L w /2 


9w{Z) = o lN ^ x (1 - e^ z ) (A13) 


2 qx (k 2 ex - q%) 
and, then, the function: 

L w /2 

~L m /2 


•“to/ “ 

(w)= J dZ* ex (Z)g w (Z) (A14) 


N 2 

M w (u) = -f- 


1 


- iq w — + 


Lii, ( ^ 

q 

l 


n 2 1 U2 _ „2 
y x "'ex j iw 


ql 2 OL - q 2 ) 

') (l — e iqxLw ) 


(A15) 


by means of that we write the new energy function 
SW(u)-,q x ): 


g0 = 


q : p s 0 (u) 


(A16) 


E 2 X - H 2 UJ 2 - i2T NR hio + q 2 p S 0 (uj)M u 
and obtain the polariton self-energy T, ex in explicit form: 

T, ex (uj : q x ) = qj S 0 (u)M w (u) (A17) 


Appendix B: UNIAXIAL CRYSTAL SLAB WITH 
IN-PLANE OPTIC C AXIS 

Let us consider an electromagnetic wave incident from 
an isotropic medium (vacuum or air) onto the left sur¬ 
face of an uniaxial multilayer. The Cartesian coordi¬ 
nate system is chosen such that the xz is the plane 
of incidence and the z axis is normal to the reflect¬ 
ing surface. The optical response of the uniaxial mul¬ 
tilayer may be characterized by four reflection ampli¬ 
tudes (R ss , R sp , R ps , R pp ) and four transmission ampli¬ 
tudes (T ss ,T S p,Tp Sl T pp ), where the second suffix refers 
to the polarization, S or P, of the reflected and transmit¬ 
ted waves, when the incident wave is in the polarization 
state of the first suffi x 25 : 26 . In order to calculate these 
amplitudes we employ the transfer-matrix method that 
connect the electromagnetic field amplitudes at both the 
interfaces of the layers. Constitutive equations for uni¬ 
axial crystal in a Cartesian framework are: 

{ (q 2 £xx ~ kf) E x + q 2 £ xy E y + (q 2 e xz + k\\k z ) E z = 0 
q £yx-k^x H - ^q £yy ^|| k z ^ E y + q £y Z E z — 0 

( q 2 Szx + k\\k z ) E x + q 2 £ zy E y + [q 2 e zz - fcfj) E z = 0 

(Bl) 


where k ex = n/L w . 
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where q = u/c and e al/ (a,u = x,y,z) are the com¬ 
ponents of the Cartesian symmetric dielectric tensor e, 
that, for uniaxial crystals with optical axis C in the re¬ 
flecting plane (xy) and forming an angle a with the x 
axis, becomes: 

/ £± + S 2 x (e|| - £x) s x s y (en - £_l) 0 

£ = s x s v (e|| - e±) ex + s 2 y (e|| - ex) 0 

\ 0 0 ex 

(B2) 

where £|| and ex are the dielectric constants, valued in 
the principal axes framework, that describe the electro¬ 
magnetic response for solicitations along the directions 
parallel and perpendicular to the optical C axis, respec¬ 
tively. s x = cos a and s y = sin a are the projections of 
C on the x and y axes. We take nonmagnetic materials 
for which the magnetic permeability value y, = 1 is used 
throughout the whole layered medium. In order to have 
nontrivial solutions of Eqs. (Bl), the determinant of the 
coefficient matrix must vanish. This requires the solution 
of the fourth degree (in k z ) characteristic equation and 
the four roots, eigenvalues of the matrix, are the wave 
vectors relative to the ordinary (±k oz ) and extraordinary 

two rays: 


k 


OZ 



(B3) 

I 


= ±1 


lk 2 e || +k\ [s 2 (e|| - ex) - £| 


£± 


(B4) 


where k 0 = (u/c) y/e±, fcy = k x = (u/c) yfe\ sin Ox. Ox 
is the angle of incidence and ei = 1. Notice that the 
expressions in Eqs. (B3) and (B4) are all independents 
of the polarization state of the incident wave. The cor¬ 
responding twelve solutions of Eq. (Bl), are Ef w e ±lkozZ 
(u = x,y,z) for the ordinary ray and E± v e lk <‘* z for the 
extraordinary ray: 


E £ } = =F Sykoz E iP = s x k 2 z (B5) 

E$ = ±s x k oz Eip = s y k 2 0 
E £ )=s y k \\ E H ] = SxkiPk^. 


Finally, each E u component (v = x, y , z) of the electric 
field is obtained as a linear combination of the corre¬ 
sponding four partial waves: 


E v = (A 0 E+ V e ik °* z + B 0 E~ v e~ ik °* z + A e E +, 




+ B e E~e ik * 


„ik\\ x 


(B6) 


The magnetic field components are then obtained by the 
Maxwell equation V x E = iufioyH. The four coeffi¬ 
cients A 0 ,B 0 ,A e and B e are obtained, for each layer of 
heterostructure, by imposing the continuity of the tan¬ 
gential components, of both the electric and magnetic 
fields, at the boundaries of each layer together to the 
following propagation relation: 


Ci 


/ E x ( Zj) \ 


E y ( z j) 

- 

H x ( Zj) 

\ Hy (Zj) y 



E x (Zj + Lj) 

E y ( z j + E j) 
{-■ ' T 


(B7) 


where Zj and Lj are the coordinate of the left interface 
and the thickness of the layer respectively. The transfer 
matrix Cj is easily building up, for each j-layer, as: 


( e 


ikioz L-i 


U(Lj) = 


g ikj oz Lj 


Jk.„ Lj 


g jez J 


(BIO) 

The transfer Matrix C of the whole n-layer stack is: 

C = C n C n -1 . C\. Finally we can write the following 

relations between the incoming and outgoing electromag¬ 
netic fields for S polarized incident wave: 


C 


cos 6\R sp 

1 + Rss 

-k\ cos 9\ (1 - R ss ) 
kxR S p 


cos 0 2 T sp 
T 

^ SS 

-k '2 cos 0 2 T SI 
k 2 T s , 


L sp 


and for P polarized incident wave: 


(BH) 


where: 


Cj = MjUj(L j )M ~ 1 


Mj = 


( E L E o X 

E oy Ky 

H+ H~ 

OX OX 

V H ty H- y 


E ex E ex \ 

E ty E ey 

H+ H~ 

ex ex 

H ty H ey J 


(B8) 


( cos Ox (1 + Rpp) \ 


/ 

cos 0 2 Tpp \ 


C 

Bps 

- 


T ps j 


kx cos OxR P s 

“ 


—k 2 cos 0 2 T ps 



\ kx (1 - Rpp ) / 


V 

k 2 Tp P ) 


(B12) 


(B9) 


where kx = q^fEx-, k 2 = q^/e^, Ox and 0 2 angles of 
incidence and transmission respectively. The refraction 
(R) and transmission (T) amplitudes are then computed 


















as the ratio of the intensities of the electric field, of the 
reflected or transmitted waves, to the incident electric 
field intensity. 
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